Zero-bias anomaly induced by the point defect in graphene 
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It is generally believed that a point defect in graphene gives rise to an impurity state at zero 
energy and causes a sharp peak in the local density of states near the defect site. We revisit the 
defect problem in graphene and find the general consensus incorrect. By both analytic and numeric 
methods, we show that the contribution to the local density of states from the impurity state 
vanishes in the thermodynamic limit. Instead, the pronounced peak of the zero-bias anomaly is a 
power-law singularity 1/\E\ from infinite resonant peaks in the low-energy regime induced by the 
defect. Our finding shows that the peak shall be viewed as a collective phenomenon rather than a 
single impurity state in previous studies. 



Graphene, a single-layer graphite composed of car- 
bon atoms arranged in two-dimensional honeycomb lat- 
tice, is recently fabricated in laboratory and attracts in- 
tense attentions from both experimental and theoreti- 
cal aspects [TJ 12 [31 H]. One of the striking features of 
graphene is its linear dispersion, allowing a relativistic de- 
scription by a pair of massless Dirac fermions in the low- 
energy regime [5J [5] . The relativistic spectrum gives rise 
to interesting phenomena such as half-integer quantum 
Hall effect [7], Klein paradox[5], edge magnetism^ ITU] 
and so on[lT][l^[L3][14]. In addition, it has been known 
for a long time that a point defect in the massless Dirac 
system induces peculiar quasi-localized state [TS]. Recent 
theoretical investigations confirm its presence as a pro- 
nounced peak in the local density of states (LDOS) (IB). 
Furthermore, the first-principles calculations seem to in- 
dicate that there exists quantized magnetic moment as- 
sociated with each defect in graphene |17l IT51 ITS] . Since a 
point defect can be created by chemisorption of hydrogen 
atoms onto graphene [TT1 |2"U1 |2"T] . these calculations sug- 
gest an alternative approach to fabricate spin qubits[22 
with relatively easy efforts, if the predicted moments sur- 
vive the quantum fluctuations not included in these stud- 
ies. 

While it is not yet settled whether the quantum fluc- 
tuations will destroy the magnetic moment near a defect, 
it is generally believed that the zero-bias anomaly in the 
LDOS corresponds to the impurity state from the point 
defect [23]. We revisit this problem with both analytic 
and numeric methods and find the general consensus on 
this issue is incorrect. We start from graphene in nan- 
otorus geometry and investigate how the system evolves 
as the size becomes infinite in the thermodynamic limit. 
At finite system size N, the defect indeed gives rise to 
a zero-energy peak in LDOS as expected |24j. However, 
its spectral weight decays as 1/lniV or l/N (depending 
on whether the nanotorus is semiconducting or metallic) 
and eventually vanishes. This conclusion, overlooked by 
previous studies, is indeed reasonable since the impurity 
state is quasi-localized and should have no contribution 



to LDOS as the system size goes infinite. 

How can we explain the zero-bias anomaly found in 
previous studies then? This paradox can be explained 
in two steps. First of all, our numerics show that the 
defect in graphene induces enormous resonant peaks in 
the LDOS at energies close to zero. Then, as the system 
size grows to infinity, these peaks crowd into zero en- 
ergy and become singular. Both numerical and analytic 
approaches give power-law singularity with weak 

logarithmic corrections. That is to say, the zero-bias 
anomaly in the LDOS is not from a single impurity state. 
Instead, it is a power-law singularity from collective reso- 
nance induced by a single defect. This is remarkable that 
the impurity state in graphene dissolve into a (weaker) 
power-law singularity as the single-particle state disap- 
pears in one-dimensional interacting electron gas. It is 
suggesting that graphene is already at criticality so that 
introduction of a point defect reshuffles the LDOS lead- 
ing to a power-law singularity rather than a well-defined 
delta-function (or broadened Lorentzian) peak. Our find- 
ing here has a significant impact on the idea of fabricating 
spin qubits by defects in graphene. Since the peak in the 
LDOS is in fact a power-law singularity from many reso- 
nant states, the quantum coherence between these states 
with realistic interactions will be tough to maintain for 
a working qubit. 

Now we walk through the details which lead to the con- 
clusions sketched in above. Because the band structure 
of graphene obtained by the first-principles calculations 
is well approximated by the nearest-neighbor hopping for 
the active ir orbitals [231 [2l>] , it is sufficient to start from 
a tight-binding Hamiltonian and add a single defect at 
the origin, 

H =~ t Yl [ ct ( r ) c ( r ') + c VMr)] + Vbct(0)c(0) , 

<r,r') 

(1) 

where t is the nearest-neighbor hopping amplitude, and 
Vq is the strength of the impurity potential. In the re- 
maining part of the calculations, we mainly focus on the 
unitary limit Vq — > oo. Though the impurity state in 
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FIG. 1: (Color online) The LDOS for a nanotorus of the 
size TV = 300 x 144 (a) with and (b) without a point defect 
evaluated at its nearest neighbor site. 
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FIG. 2: (Color online) The height of an individual resonant 
peak in the LDOS close to the zero energy versus the system 
size TV. The green dashed line is the fit to the 1/lnN scaling. 



the unitary limit has been solved analytically in a recent 
paper [16^ it is insightful to rederive it with cares so that 
the evolution of the impurity state with the system size 
TV is clarified. We follow the definitions used in Ref. [16] 
and reformulate the problem in more compact notations. 

Consider a zigzag carbon nanotube with TV unit cells 
around the transverse direction. The total number of 
localized states at the left (right) of the defect is denoted 
as Nl(Nr). Construct a (Nl + iV^)-dimensional vector 
to represent the solution, A = [a^,(l + e %km ')a k J, 
Meanwhile, we can also introduce (TV — 1) vectors in the 
same space Bj = [e lfemJ , e ikn >' 3 ], where j — 1,2, TV — 1. 
The boundary conditions |16j then take the simple form 
A Bj = 0, i.e. we are looking for all linearly independent 
vectors A which is orthogonal to the subspace expanded 
by the TV — 1 vectors Bj . 

For semiconducting nanotubes (TV ^ 3n), the quan- 
tized momentum does not cut through the Dirac points 
at k — 27r/3,47r/3. Thus, counting all localized states 
within 2tt/3 < k m < 4tt/3 gives N L = [(TV+l)/3], where 
[x] denotes the Gauss symbol. Similarly, it is straight- 
forward to obtain Nr = [(2TV + l)/3]. Since the com- 
bination of k m and k m i exhausts all quantized momenta 
in the Brillouin zone, TVl + Nr = N. The boundary 
conditions nail down the symmetric solution, 

A = [1, 1, 1]. (2) 

It is easy to verify that A B 3 = ~£n=i e l2 ™ j / N = 
NSj.o = 0, satisfying all orthogonal criteria. Further- 
more, because Bj expand a subspace of dimension N — 1 
in the TV dimensional space, only one solution is allowed. 
Normalization of the above solution involves a summa- 
tion over all localized states, ^ m [l — 2cos(fc m /2)] ~ 
J dk/k ~ In TV. As a result, the spectral weight at E = 
scales as 1/ In TV in the thermodynamic limit. 

For metallic nanotubes (N = 3n), there are two addi- 
tional extended states at the Dirac points with the nodal 



structure. These extended states can hybridize with the 
localized states and should be included on both sides. 
Thus, Nl = (N/3) + 1 and N R = (2JV/3) + 1. The di- 
mension of the vector space is larger Nl + Nr — N + 2 
here. Let us count the number of independent solutions 
first. Subtracting the dimension of the vector space TV +2 
by the dimension TV — 1 expanded by Bj , we expect three 
independent solutions. Let us start with the simple ones 
which do not involve the localized states, 

A m = [1,0, 0,...,0, -1,0, 0,...,0 ], 

A D2 = [0,1, 0,...,0, 0,-1, 0,...,0 ]. (3) 

For clarity, the components corresponding to the Dirac 
points are in italics while those for localized states are in 
roman. The orthogonal criteria are trivially satisfied. In 
fact, these two solutions are just the extended states at 
the Dirac points. The third solution is 



Although the solution is similar to that in Eq. ([2|, its 
hybridization with the extended state at the Dirac points 
makes it an extended state as well. Thus, despite the 
superficial resemblance, its contribution to the spectral 
weight at E — is dramatically different and scales down 
much faster as 1/TV. 

The above analytic solutions exclude the impurity 
state as the cause for the zero-bias anomaly of the LDOS 
in graphene. We also perform numerical diagonaliza- 
tion to confirm the different trends for semiconducting 
and metallic cases. Note that , in addition to the impu- 
rity state, the defect in graphene also generates resonant 
states at energies close to E = inaccessible by ana- 
lytic approach. We employ numerical diagonalization of 
the carbon nanotorus (by jointing the two ends of the 
zigzag carbon nanotube together) to address this issue. 
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For comparison, we computed the LDOS for a metal- 
lic nanotorus with and without a point defect in Fig. [T| 
The change of the LDOS right at zero energy caused by 
the defect is almost invisible. But, the induced resonant 
peaks close to the zero energy are enormous and largely 
change the profile of the LDOS. It is reasonable to ex- 
pect that the resonant peaks are far more important than 
the impurity state at zero energy in the thermodynamic 
limit. 

Furthermore, our numerics also reveal that the height 
of each individual resonant peak at finite (but small) en- 
ergy scales as 1/lniV with the system size as shown in 
Fig. [2] This 1/ In N dependence of the resonant peaks 
bridge a gap between the semiconducting and the metal- 
lic nanotori in thermodynamic limit. Note that the impu- 
rity state by itself in the semiconducting or in the metallic 
nanotorus gives 1/ In N or 1 /N contribution to the zero- 
energy LDOS. The difference in LDOS due to the impu- 
rity state in the two classes is huge in a large but finite 
system. The merger of the 1/ In iV-resonant peaks re- 
solves this discrepancy between the semiconducting and 
metallic systems and establishes the smooth thermody- 
namic limit to the two-dimensional graphene. 

The remaining task is to understand how these res- 
onant peaks crowd into the zero-energy regime as the 
system size goes to infinity. The numerical results are 
shown in Fig. [3] In the inset, we compare the LDOS 
near the defect and in the bulk. As in previous stud- 
ies, the most noticeable difference is the sharp peak at 
E = 0. However, one may also observe that there ex- 
ists large spectral weight transfer from the higher-energy 
regime to the lower energy. This is an indirect hint that 
the sharp peak at zero energy may not attribute to a 
single impurity state. To illustrate this point, the LDOS 
shown in the log-log plot reveals the 1/\E\ power-law de- 
pendence. The deviation from the power-law singularity 
at extremely small energy where it is rounded off comes 
from the Lorentzian energy broadening factor introduced 
in the numerical calculations. 

To further strengthen our numerical findings, we solve 
the Dyson equation for graphene in the presence of a 
point defect. The change of the LDOS can be expressed 
in terms of the non-interacting retarded Green's func- 
tions, Gq (k,E) = 1/(E ± \h(k) \ + ir)), representing the 
propagation of particles and holes with the dispersion 
h(k) = t + 2tcos(k x /2)e i ^ k y/ 2 from the tight-binding 
model. Following standard Green's function techniques, 
the change of the LDOS at the first-nearest-neighbor sites 
(B) of the defect (located at the A site) is 



Ap(E) = --Im 

7T 



VqGabGba 
1 - V Gaa 



(5) 



where Gaa' with A, A' = A, B are the Green's functions 
between sites A and A'. In the unitary limit Vo — > 00, the 
change of the LDOS no longer depends on the strength 
of the impurity potential as expected. 
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FIG. 3: (Color online) The LDOS at the nearest-neighbor 
site of the defect for a nanotorus of the size N = 564 x 564 
with a broadening factor 0.004 (red circles) and in the two di- 
mensional limit with a broadening factor 10 _J (blue squares). 
The green dashed line of 1/\E\ is the guide to the eyes. In the 
inset, the LDOS at the nearest-neighbor site (shown in red 
solid line) is compared to that deep inside the bulk (shown in 
blue dashed line). Since the peak at zero energy is high, it is 
chopped off by the double-arrow dashed line for clarity. 



For graphene, it is straightforward to show that sym- 
metries between these Green's functions lead to the rela- 
tions, Gaa = G B b and Gab = G B a, 



Gaa(E) = |0[G o -(k,B) + G+(k ) £)] ) 
d 2 k h(k) 



(0) 



Gab(E) 



8tt 2 \h(k)\ 



[G -(k,£)-G+(k, £)].(?) 



These Green's functions can be solved numerically to 
compute the LDOS. The results are shown as blue 
squares in Fig. [3] It delivers the same power-law singu- 
larity in the LDOS. It is also interesting to notice that, 
though with the same exponent, the absolute values of 
the LDOS for the large nanotorus with N — 564 x 564 
and the two-dimensional graphene are different - another 
hint for the collective phenomena rather than a single im- 
purity state for the zero-bias anomaly. 

The power-law singularity can be derived from the lin- 
ear dispersion near the Dirac cones. Carrying out the an- 
gular part of the integral near the Dirac points with linear 
dispersion, the Green's functions are approximately 



Gaa(E) ~ Eln(E 2 /A 2 )+i\E\. 
Gab(E) ~ A 2 - E 2 ln(E 2 /A 2 ) 



(8) 

isign(E)E 2 , (9) 



where A is a momentum cutoff introduced for lineariz- 
ing the spectrum. Since Gab(E) appears in the nu- 
merator, it can be treated as a constant in low-energy 
limit. Therefore, the change of the LDOS is Ap(E) ~ 
Im[GAA]/(Re[GAyi] 2 + Li^G^a] 2 ) an d the singularity 
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near E = emerges, 
Ap(E) 



|£|(ln|E|r 



(10) 



Because the logarithmic correction is very weak (beyond 
the resolution of our numerical results), the singularity 
is essentially a power-law l/\E\ and agrees with our pre- 
vious numerical findings. 

The above calculations are readily generalized to other 
systems with gapless dispersion h(k) ~ |k|". Follow- 
ing the same steps, it is straightforward to show that 
the imaginary of the Green's function follows the energy 
dependence, lm[G AA \ ~ |^|( 2 /«)-i. I gnoring the loga- 
rithmic correction, the real part of the Green's function 
shares the same energy dependence. Thus, the change of 
the LDOS exhibits the power-law Ap(E) oc \E\ 1 ~ {2 / a \ 
For double-layer graphene, the dispersion is quadratic 
with a = 2, giving a constant Ap(E) at zero energy. Our 
calculations thus predict that the power-law singularity 
is absent for the bilayer graphene in sharp contrast to the 
1/\E\ singularity for the single-layer. 

So far, no correlation effect is taken into account. 
Though it is generally believed that the mutual inter- 
actions between electrons are likely to be irrelevant in 
graphene, there are evidences that the correlation ef- 
fects may be significant near defects or open bound- 
aries. Further in-depth investigations are necessary to 
explore how the power-law singularity evolves with the 
inclusion of electronic correlations. However, it is worth 
mentioning how the exponent changes in the analogous 
one-dimensional interacting system. Upon the inclusion 
of interaction, the delta-function for the quasi-particle 
changes into 1/|-E| power- law singularity. Increasing 
the strength of the interaction suppresses the density of 
states near the zero energy and gradually changes the 
exponent from the negative 7 — —1 to positive values, 
exhibiting the so-called "pseudo-gap" behavior. We may 
expect similar trend happens for defect in graphene but 
it needs further investigations. 

In conclusions, we revisit the defect problem in 
graphene and find the sharp peak in the LDOS is not 
due to the impurity state at zero energy. Instead, it is a 
power-law singularity l/\E\ from infinite resonant peaks 
in the low-energy regime induced by the defect. Our find- 
ing here has a significant impact on the idea of fabricating 
spin qubits by defects in graphene. Since the peak in the 
LDOS is in fact a power-law singularity from many reso- 
nant states, the quantum coherence between these states 
with realistic interactions will be tough to maintain for 
a working qubit. 
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